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ABSTRACT 

Detection of the radiation emitted from the first galaxies at z > 10 wiU be made pos- 
sible in the next decade, with the launch of the James Webb Space Telescope (JWST). 
We carry out cosmological radiation hydrodynamics simulations of Population III (Pop 
III) starbursts in a 10* M© dwarf galaxy at z ~ 12.5. For different star formation ef- 
ficiencies and stellar initial mass functions (IMFs), we calculate the luminosities and 
equivalent widths (EWs) of the recombination lines Ha, Lya, and He ii A1640, under 
the simple assumption that the stellar population does not evolve over the first 3 
Myr of the starburst. Although only < 40 percent of the gas in the central 100 pc of 
the galaxy is photoionized, we find that photoheating by massive stars causes a strong 
dynamical response, which results in a weak correlation between luminosity emitted 
in hydrogen recombination lines and the total mass in stars. However, owing to the 
low escape fraction of He ii-ionizing photons, the luminosity emitted in He ii A1640 is 
much more strongly correlated with the total stellar mass. The ratio of the luminosity 
in He ii A1640 to that in Lya or Ha is a complex function of the density field and the 
star formation rate, but is found to be a good indicator of the IMF in many cases. 
The ratio of observable fluxes is FAi64o/FHa ^ 1 for clusters of 100 Pop III stars 
and FAi64o/FHa 0.1 for clusters of 25 Mq Pop III stars. The EW of the He ii A1640 
emission line is the most reliable IMF indicator, its value varying between ^ 20 and 
~ 200 A for a massive and very massive Pop III IMF, respectively. Even the bright, 
initial stages of Pop III starbursts in the first dwarf galaxies will likely not be directly 
detectable by the JWST, except in cases where the flux from these galaxies is strongly 
magnified through gravtitational lensing. Instead, the JWST may discover already 
more massive, and hence more chemically evolved, galaxies in which primordial star 
formation has largely ceased, or is contaminated with more normal, Pop I/II, star 
formation. 

Key words: cosmology: theory - early Universe - galaxies: formation - H ii regions 
- ultraviolet: galaxies. 



1 INTRODUCTION 

The epoch of the first galaxies marked a fundamental tran- 
sition in the Universe, ending the Cosmic Dark Ages, begin- 
ning the process of reionization, and witnessing the rapid 
proliferation of star formation and black hole growth (e.g. 
Barkana & Loeb 2001; Bromm & Larson 2004). While the 
theory of primordial star formation and early galaxy for- 
mation has rapidly developed (see Ciardi & Ferrara 2005; 
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Glover 2005; Barkana & Loeb 2007), observations of the 
first galaxies at redshifts z ^ 10 have so far been out of 
reach (but see Stark et al. 2007). In the coming decade, the 
James Webb Space Telescope (JWST) promises to provide 
direct observations of this critical period in cosmic history, 
allowing to place new constraints on the stellar initial mass 
function (IMF) at high redshift, on the luminosity function 
of the first galaxies, and on the progress of the early stages 
of reionization (e.g. Barton et al. 2004; Gardner et al. 2006; 
Windhorst et al. 2006; Ricotti et al. 2008; Haiman 2008). 
The IMF of the stellar populations which form in the 
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first galaxies is of central importance in determining their 
properties and impact on early cosmic evolution. The cur- 
rent theoretical consensus posits that the first stars, which 
formed in isolation in dark matter (DM) minihaloes, likely 
had masses of the order of 100 Mq (Bromm et al. 1999, 
2002; Abel et al. 2002; Yoshida et al. 2006, 2008; McKee & 
Tan 2008). In the first galaxies, which form in DM haloes 
with masses ^ 10^ Mq (Greif et al. 2008; Wise et al. 2008), 
there is no such theoretical consensus on the IMF, as the 
initial conditions of the star-forming gas are uncertain (see 
Jappsen et al. 2009a, b) . A large fraction of these first galax- 
ies are likely to already host Population II (Pop II) star 
formation, due to previous metal enrichment (Johnson et 
al. 2008; see also Trenti & Stiavelh 2009; Clark et al. 2008). 
However, clusters of primordial stars likely form in some 
fraction of the first galaxies, owing to either inhomogeneities 
in the Lyman- Werner (LW) background which can suppress 
star formation where galaxies are strongly clustered (Dijk- 
stra et al. 2008; see also Ahn et al. 2008) or to the direct 
collapse of the first stars to black holes, thereby locking up 
the metals produced in their cores (e.g. Heger et al. 2003). 

Due to the hard spectra of massive metal-free stars, 
strong nebular emission in helium recombination lines has 
been suggested as an observable indicator of a population of 
such stars (e.g. Bromm et al. 2001; Oh et al. 2001; Tumlin- 
son et al. 2001; Schaerer 2002). In particular, a high ratio 
of the luminosity emitted in He ll A1640 to that emitted in 
Lya or Ha may be a signature of a galaxy hosting massive 
Pop HI star formation, and has already served as the basis 
for searches for such galaxies (e.g. Dawson et al. 2004; Na- 
gao et al. 2005, 2008). In addition, high equivalent widths 
(EWs) of Lya and He ll A1640 are expected to character- 
ize galaxies undergoing a Pop HI starburst (e.g. Schaerer 
2003). While no definitive detections of such galaxies have 
been achieved to date, observations of galaxies at 3 ^ 2 ^ 
6.5 with large Lya EW and strong He ll A1640 emission may 
indicate that some galaxies host Pop HI star formation even 
at such relatively low redshift (see Jimenez & Haiman 2006; 
Dijkstra & Wyithe 2007). 

Previous analytical calculations of the recombination 
radiation expected from the first galaxies have been car- 
ried out under a number of idealized assumptions, namely 
of a static, uniform density field, and of the formation of a 
static Stromgren sphere. Taking a complementary approach, 
we study here the properties of the recombination radia- 
tion emitted by the first galaxies with a focus on how the 
dynamical evolution of the galaxy affects the properties of 
this radiation. We present high-resolution cosmological radi- 
ation hydrodynamics simulations of the production of nebu- 
lar emission from a cluster of primordial stars formed within 
the first galaxies. We resolve the H ll and He ill regions gen- 
erated by the stellar cluster, thereby arriving at improved 
predictions for the emission properties of the first galaxies, 
which will be tested by the JWST. 

Our paper is organized as follows. In Section 2, we de- 
scribe our simulations and the methods used in their anal- 
ysis; in Section 3 our results are presented, along with the 
implications for both IMF and star formation rate (SFR) in- 
dicators; in the final Section 4, we summarize and give our 
conclusions. 




Figure 1. The projected gas temperature ai z = 12.7, just before 
the stellar cluster turns on. Shown at left is the entire 1.4 Mpc 
(comoving) simulation box. The region of highest resolution is in 
the center of our multi-grid simulation box; magnified at right 
is the central ~ 10 kpc (physical), where the ~ 10* Mq dark 
matter halo hosting the stellar cluster is located. Note that the 
gas is heated to > 10^ K when passing through the virialization 
shock. 

2 METHODOLOGY 

We carry out radiation hydrodynamics simulations which 
track the impact of the radiation from Pop HI stellar clusters 
forming promptly within a dwarf galaxy at z ~ 12.5. In 
this section we describe the simulations and the calculations 
carried out in analyzing their outcomes. 

2.1 The simulations 

As with previous work, for our three-dimensional numer- 
ical simulations we employ the parallel version of GAD- 
GET (version 1), which includes a tree (hierarchical) gravity 
solver combined with the smoothed particle hydrodynamics 
(SPH) method for tracking the evolution of gas (Springel et 
al. 2001; Springel & Hernquist 2002). Along with H2, H+, H, 
H~, H^, e~. He, He"*", and He"*""*", we have included the five 
deuterium species D, D"*", D~, HD and HD^, using the same 
chemical network as in Johnson & Bromm (2006, 2007). 

For our simulation of the assembly of a dwarf galaxy 
at z ~ 12.5, we have employed multi-grid initial conditions 
which offer higher resolution in the region where the galaxy 
forms (e.g. Kawata & Gibson 2003). We initialize the simu- 
lation according to the ACDM power spectrum at z — 100, 
adopting the cosmological parameters fim = 1 — ~ 0.3, 
Sis = 0.045, h = 0.7, and as = 0.9, close to the values mea- 
sured by the Wilkinson Microwave Anisotropy Probe in its 
first year (Spergel et al. 2003). Here we use a periodic box 
with a comoving size L = 1 Mpc for the parent grid. 
Our simulations use A^dm ~ A^sph = 1.05 x 10® particles for 
DM and gas, where the SPH particle mass is mspH ~ 120 
Mq in the region with the highest resolution. For further de- 
tails on the technique employed to generate our multi-grid 
initial conditions, see Greif et al. (2008). The maximum gas 
density that we resolve is Urcs ~ 10'^ cm"'', while gas at 
higher densities is accreted onto sink particles, as described 
in Johnson et al. (2007) . We have also included the effect of 
a LW background radiation field, at a level of Jlw = 0.04 
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X 10 erg s ^ cm ^ Hz ^ sr ^ , just as in Johnson et al. 
(2008). 

To capture the effects of the ionizing radiation emitted 
by a single Pop III stellar cluster within the dwaxf galaxy, 
we approximate the cluster as a point source located at the 
center of the most massive DM halo in our simulation box at 
z = 12.7. The projected gas temperature in the entire sim- 
ulation box at this redshift is shown in Fig. 1 (left panel); 
the right panel shows the temperature in the region of the 
host halo at the center of the box, which has a virial mass 
of 9 X lO'^ Mq, characteristic of the first galaxies. At each 
timestep, we find the boundaries of both the H ii and He iii 
regions generated by the stellar cluster using a ray-tracing 
technique that improves our earlier implementation (John- 
son et al. 2007). 

The procedure used to calculate the Stromgren sphere 
around the stellar cluster for a given time-step At is similar 
to the ray-tracing scheme used in Johnson et al. (2007). 
We create a spherical grid centered at the location of the 
cluster, consisting of ~ 1.2 x 10'' rays and 1000 linearly 
spaced radial bins. We resolve the central kiloparsec around 
the source, roughly the virial radius of the host halo, with 
250 radial bins, while the remaining 750 bins are linearly 
spaced out to ~ 20 kpc. 

In a single, parallelized loop, the Cartesian coordinates 
of all particles are converted to spherical coordinates, such 
that their density and chemical abundances may be mapped 
to the bins corresponding to their radius, azimuth and 
zenith, denoted by r, 9 and (j), respectively. To avoid miss- 
ing dense clumps, particles contribute to bins independent 
of distance, but proportional to their density squared. Once 
this preliminary step is complete, it is straightforward to 
solve the ionization front equation along each ray: 

2dri Qion Z"'^' 2, 

n„ri — - = — QB / Hen+r dr , (1) 

dt 47r Jo 

where n denotes the position of the ionization front, Qion 
the number of ionizing photons emitted per second by the 
stellar cluster, qb the case B recombination coefficient, and 
Un, rie and n+ the number densities of nonionized particles, 
electrons and ionized particles, respectively. The numbers 
of H I-, He I- and He ll-ionizing photons arc Qion = A^*9ion, 
where A'', is the number of stars in the cluster (here we 
assume that all have the same mass) and gion is the number 
of ionizing photons emitted by a single star, given by 

where a is the Stefan-Boltzmann constant, Umin denotes 
the minimum frequency corresponding to the ionization 
thresholds of H I, He I and He ii, and we assume that 
massive Pop III stars emit a blackbody spectrum Bi, (in 
erg s~' cm~^ Hz"'^ sr~^) with an effective temperature Tgg 
and a luminosity L, (e.g. Schaerer 2002). 

To obtain a discretization of the ionization front equa- 
tion, we replace the integral on the right-hand side of equa- 
tion (1) by a discrete sum: 
/■n 

/ UeU+r^dr = ^ rie,in+,irf Ar , (3) 
Jo 

where Ar is the radial extent of the individual bins. Simi- 
larly, the left-hand side of equation (1), which models the 



propagation of the ionization front into neutral gas, is dis- 
cretized by 

"nr-f^ = ■^'^nr,,iriAri , (4) 

i 

where At is the current time-step and the summation is over 
radial bins starting with the bin lying immediately outside 
of ri,oid, the position of the ionization front at the end of 
the previous time-step, and ending with the bin lying at the 
new position of the ionization front. We perform the above 
steps separately for the H ii and He iii regions, since they 
require distinct heating and ionization rates. 

We have chosen the size of the bins that are used in our 
ray-tracing routine to roughly match the volume of gas rep- 
resented by a single SPH particle within the ^ 1 kpc virial 
radius of the halo hosting the stellar cluster, such that the 
boundaries of the photoionized regions are maximally re- 
solved while also reliably conserving ionizing photons. How- 
ever, in some cases it may occur that the mass contained 
in a bin is smaller than that of the SPH particle contained 
within it, such that ionizing the entire SPH particle involves 
ionizing more gas than is contained in the bin. In turn, this 
can lead to an overestimate of the number of recombina- 
tions. While this effect is minor in our simulations, in the 
calculations presented below we enforce that the total num- 
ber of recombinations does not exceed the total number of 
ionizing photons available. 

We carry out four simulations, each with a different 
combination of IMF and total cluster mass. For the IMF, 
we assume for simplicity, and in light of the still complete 
uncertainty regarding its detailed shape, that the cluster 
consists either entirely of 25 Mq or 100 Mq Pop III stars. 
These choices are meant to bracket the expected charac- 
teristic mass for Pop III stars formed in the first galaxies, 
which depending on the cooling properties of the gas may 
be Pop III.2 stars with masses of order 10 Mq or, possibly. 
Pop III.l stars with masses perhaps an order of magnitude 
higher (Johnson & Bromm 2006; Greif et al. 2008; McKee 
& Tan 2008; but see Jappsen et al. 2009a). For each of these 
IMFs, we vary the total stellar mass in the cluster, choosing 
either 2.5 x 10^ Mg or 2.5 x 10** Mq for the total mass in 
stars. These choices correspond to ^ 1 and ~ 10 percent, re- 
spectively, of the cold gas available for star formation within 
the central few parsecs of such a primordial dwarf galaxy 
(see Wise et al. 2008; Regan & Haehneh 2009). We calcu- 
late the ionizing flux from each of these clusters, assuming 
blackbody stellar spectra at 7 x 10'' and 10^ K and bolo- 
metric luminosities of 6 x 10'' and 10^ L©, for the 25 and 
100 Mq stars, respectively, appropriate for metal-free stars 
on the main sequence (Marigo et al. 2001). 

For simplicity, we have chosen to keep the input stellar 
spectra constant in time over the course of the simulations. 
Accordingly, we run the simulations only for 3 Myr, which 
is roughly the hydrogen-burning timescale of a 100 Mq pri- 
mordial star, and about half that of a 25 Mq primordial 
star. We note that while the H l-ionizing flux from 100 Mq 
primordial stars is roughly constant over this timescale, the 
He ii-ionizing flux decreases by a factor of ~ 2 by a stellar 
age of 2 Myr, and by a much larger factor near the end of 
hydrogen-burning as the star evolves to the red (e.g. Marigo 
et al. 2001; Schaerer 2002). We note, however, that stel- 
lar models accounting for the effects of rotation yield less 
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precipitous drops in the emitted ionizing fiux with time, as 
fast rotation, especially of low-metallicity stars, can keep 
the stars on bluer evolutionary tracks (e.g. Yoon & Langer 
2005; Woosley & Heger 2006; Vazquez et al. 2007); indeed. 
Pop III stars may have been fast rotators (see Chiappini et 
al. 2008). Nonetheless, the results that we derive pertaining 
to He II recombination emission from clusters of 100 
stars may be, strictly speaking, only reliable for stellar ages 
< 2 Myr. An in-depth study of the impact that stellar evo- 
lution has on the emission properties of primordial galaxies 
is given in Schaerer (2002); in the present work, we take a 
complementary approach and instead focus on how the emis- 
sion properties are affected by the hydrodynamic evolution 
of the gas in the first galaxies. 

We make the related simplifying assumption that the 
stellar cluster forms instantaneously. This is valid if the 
timescale for star formation tsF is much shorter than the 
lifetime of the stars that we consider, or ^sf ^ 3 Myr. If we 
assume that stars form on the order of the free-fall time iff, 
and take it that the star cluster forms within the central ~ 
1 pc of the halo (see e.g. Wise et al.2008; Regan & Haehnelt 
2009), then we find ts¥ ~ 5 x 10® yr, for which our assump- 
tion is marginally valid. We note that more work is needed 
to accurately determine the star formation timescale in the 
first galaxies, as the works cited here neglect, in particular, 
the important effect of molecular cooling on the evolution 
of the primordial gas. 

2.2 Deriving the observational signature 

The simulations described above allow us to calculate the lu- 
minosities and equivalent widths of the recombination lines 
emitted from high-redshift dwarf galaxies during a primor- 
dial starburst. A related quantity we obtain is the escape 
fraction of ionizing photons from such a galaxy. Here we 
describe each of these calculations. 



2.2.1 Escape fraction of ionizing photons 

Photons which escape the host halo from which they are 
emitted proceed to reionize the intergalactic medium (IGM), 
where densities are generally very low, yielding long recom- 
bination times. Ionizing photons which do not escape the 
host halo are, however, available to ionize dense gas which 
recombines quickly, leading to appreciable emission in re- 
combination lines. Therefore, the luminosity of a galaxy in 
recombination radiation is intimately related to the escape 
fraction of ionizing photons. The escape fraction of ionizing 
photons from the halo hosting the stellar cluster is given by 
subtracting the number of recombinations Qrcc per second 
within the virial radius from the total number of ionizing 
photons emitted by the cluster: 

/esc — (Qion Qrec) (Qion) , (5) 

again with Qion ~ Af*9ion. This equation is valid under the 
assumption that within the host halo the number of ioniz- 
ing photons which fail to escape is balanced by the number 
of recombinations within the halo. This is a reasonable as- 
sumption, since the number of atoms which become ionized 
within the host halo is far less than the total number of 
recombinations that occur in the halo, the ionization of pre- 
viously neutral gas being the only other sink for ionizing 



photons within the halo. The number of recombinations is 
given as 

n 2 



as — 
Pi 



/e/H 



(6) 



where q;b is the case B recombination coefficient for hydro- 
gen, TTiH the mass of a hydrogen atom, while rrii, /if and pi 
are the total mass, mean molecular weight, and mass density 
of the ith SPH particle, respectively. For each SPH particle, 
/hii and /e denote the fraction of nuclei in H ii and the 
fraction of free electrons, respectively. 

Here the summation is over all SPH particles within the 
virial radius of the host halo, or within a physical distance 
of ~ 1 kpc from the central stellar cluster. We calculate 
the escape fractions of both H l-ionizing and He ll-ionizing 
photons. These quantities are generally not equal, and they 
each contribute to determining the radiative signature of the 
initial starbursts in the first galaxies. 



2.2.2 Luminosity in recombination lines 

For each of our simulations, we compute the luminosity emit- 
ted from photoionizcd regions in each of three recombina- 
tion lines: Ha, Lya, and He II A1640. These luminosities arc 
calculated by again summing up the contributions from all 
SPH particles within the virial radius, where virtually all of 
the recombination line luminosity emerges, as follows: 
n 2 

Pi 



Pi 



/iimH 



E 



mi 
Pi 



/iimn 



/e/mi 



/e/mi 



.7A1640- 



Pi 



/e/Helll , 



(7) 



(8) 



(9) 



where the j are the temperature-dependent emission coeffi- 
cients for the lines (Osterbrock & Ferland 2006), and /Heiii is 
the fraction of helium nuclei in He III for each SPH particle. 

Given the luminosity in a recombination line over an 
area of the sky, we may compute the flux in that line, as ob- 
served at z = with a spectral resolution R = A/ A A, where 
A is the wavelength at which the emission line is observed 
(e.g. Oh et al. 2001). While Lya photons are scattered out 
of the line of sight in the IGM prior to reionization (Loeb 
& Rybicki 1999), a process which we do treat in the present 
calculations. Ha and He ll A1640 photons will not suffer such 
severe attenuation. Assuming that the line is unresolved, the 
monochromatic flux in Ha, for example, is 



/ho 



^HaAHa(l + z)R 

A-KcDl{z) 



2°"J^(io4Tgs-0(^) (i^) 



(10) 



where luct is the luminosity in Ha along the line of sight 
through the emitting galaxy, D-l{z) is the luminosity dis- 
tance at redshift z (~ 10^ Gpc at z — 10), and Ahc« is the 
rest frame wavelength of the line, 656.3 nm. If the galaxy 
is spatially unresolved, appearing as a point source, we may 
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simply substitute Lhh for fea in equation (10), to compute 
the total flux from the galaxy. In terms of total (integrated) 
line flux, we have the equivalent expression: 



7.5 kpc 



7.5 kpc 



1 kpc 



10 erg s 



1040ergs- 



l + z 
10 



2.2.3 Recombination line equivalent widths 

Another observable quantity obtained from our simulations 
is the rest-frame equivalent width (EW) of recombination 
lines. We calculate the EWs of the three recombination lines 
considered, following Schaerer (2002): 

Lua 



W 



L\,ncb + Lx,* 
_ ^A1640 



(11) 



(12) 



(13) 



where the monochromatic continuum luminosity, evaluated 
at the wavelength of the line, is the sum of the nebular 
emission L^.nob and the stellar emission La,*- The nebular 
continuum luminosity is given by 

i'A.nob = To "^"""Qrcc , (14) 

Ob 

where A is the wavelength in question, and Qroc is again 
the total number of recombinations per second in the halo. 
The continuous emission coefficient 7tot accounts for free- 
free, free-bound, and two-photon continuum emission, as 
described in Schaerer (2002). The stellar continuum lumi- 
nosity is calculated assuming a blackbody stellar spectrum 
and is given by 



(15) 



AS exp(/ic/AfcBreff) - 1 

where A'', is the number of stars in the cluster, TcB is the 
effective surface temperature of a star, and 7?, is the stellar 
radius. 



3 RESULTS AND IMPLICATIONS 

We next discuss the observable characteristics of primordial 
dwarf galaxies. In particular, we evaluate the utility of indi- 
cators for the SFR and the stellar IMF in such galaxies. 



3.1 Evolution of gas inside the galaxy 

With the ignition of a stellar cluster at the center of the 
host halo, the gas surrounding the cluster is photoheated, 
raising its pressure and leading to its outward expansion. In 
turn, the overall recombination rate in the host halo drops, 
allowing the expansion of the H ll region to continue for 
a constant rate of ionizing photon production. Fig. 2 shows 
the growth of the H ll region and the concomitant expansion 
of the gas in the center of the host halo for the more massive 
100 Mq cluster. The H ll region breaks out of the host halo 
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Figure 2. The density- weighted temperature (left column), 
density-weighted H II fraction (middle column), and number den- 
sity (right column), each averaged along the line of sight, of the 
gas surrounding the more massive 100 Mq star cluster, shown at 
three different times from the prompt formation of the cluster: 
500,000 yr (top row), 1 Myr (middle row), and 3 Myr (bottom 
row). The H II region grows as the density of the gas in the center 
of the host halo gradually drops in response to the intense pho- 
toheating. Note the different length scale of each column, given 
at the top in physical units; the density is shown only within the 
central region of the host halo. 



within the first ^ 1 Myr, and after 3 Myr it extends to ^ 
7 physical kpc, only slightly larger than the size of the H ll 
region created by a single massive Pop III star in a minihalo 
(e.g. Alvarez et al. 2006). 

While the gas surrounding the formation sites of the 
first stars in minihaloes is easily photoevacuated by a single 
massive Pop III star (e.g. Kitayama et al. 2004; Whalen 
et al. 2004), the deeper gravitational potential well of the 
DM haloes hosting the first galaxies allows for the retention 
of gas even under intense photoheating; indeed, this is one 
criterion used to define the first galaxies (e.g. Read 2006; 
Johnson et al. 2008; Greif et al. 2008). 

As shown in Fig. 3, a substantial portion of the gas 
in the galaxy, even within ~ 100 pc of the stellar clus- 
ter, remains neutral after 3 Myr. This gas is shielded from 
the ionizing radiation, causing the ionization front (I-front) 
to propagate outward anisotropically in the inhomogeneous 
cosmological density field (see also e.g. Shapiro et al. 2004; 
Abel et al. 2006; Alvarez et al. 2006). Even for the case of 
the highest ionizing flux, the fraction of ionized gas within 
the central ~ 100 pc is ^ 0.4, leaving the majority of the 
high density gas neutral. While the photodissociating radi- 
ation from the initial stellar cluster will slow the collapse 
of this primordial gas (e.g. Susa & Umemura 2006; Ahn & 
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100 1000 
r [pc] 

Figure 3. The mass fraction of enclosed gas wliicli is photoion- 
ized, /ion; ^ a- function of the distance from the central star 
cluster, after 3 Myr of photoheating. Each line corresponds to a 
different choice of IMF and total mass in stars, as labeled. The 
ionized fraction begins to drop at > 300 pc mostly due to two 
minihaloes at ~ 400 and ~ 500 pc, shown in Fig. 2 (right panels), 
which remain shielded from the ionizing radiation and are thus 
largely neutral. Within ~ 100 pc, there is a large fraction of gas 
which remains un- ionized at densities > 100 cm~^. This gas will 
likely collapse to form stars, despite the strong radiative feedback 
from the central stellar cluster. 



Shapiro 2007; Whalen et al. 2008), some fraction of it will 
likely be converted into stars once the most massive stars in 
the cluster have died out. Indeed, the shocks engendered by 
the supernovae that mark the end of their lives may expedite 
the collapse of the gas (e.g. Mackey et al. 2003; Salvaterra 
et al. 2004; Machida et al. 2005; Greif et al. 2007; Sakuma 
& Susa 2009). The incomplete ionization of the central gas 
confirms that the masses that we have chosen for the clus- 
ters are not overly large, as there is still some neutral gas 
available for subsequent star formation regardless of the ra- 
diative feedback. 

The gas that is photoionized, however, is gradually ex- 
pelled from the center of the halo, and after 3 Myr of photo- 
heating the density of the photoionized gas within ~ 20 pc 
of the cluster drops to ^ 10 cm~^ for the more massive clus- 
ter of 100 M© stars shown in Fig. 2. For our other choices of 
IMF and total cluster mass the dynamical response is less 
dramatic, as the ionizing flux is weaker; for example, after 
3 Myr the density of the central photoionized gas is ^ 50 
cm~^ for the less massive cluster of 25 M0 stars. The vary- 
ing degree to which photoheating dynamically impacts the 
host halo leads to important differences in the properties of 
the emitted radiation. 

Although the limited resolution of our simulations al- 
lows only to track the expansion of the H 11 region from an 
initial physical size of ^ 10 pc, we expect that after roughly 
a sound- crossing time of the central unresolved ~ 10 pc, or 



after the first few 10^ yr, the evolution of the H 11 region is 
reliably resolved. It should thus be noted that the breakout 
of the H II region may be delayed by of order this timescale 
compared to our simulations. We note that in the Milky Way 
the expansion of the photoheated gas in an H 11 region may 
be slowed due to turbulent pressure confinement (e.g. Xie et 
al. 1996; Mac Low et al. 2007), resulting in ultra-compact 
H II regions persisting for ^ 10^ yr (Wood & Churchwell 
1989), much longer than the sound-crossing timescale for 
such regions. As turbulence begins to play an important role 
in the formation of the first galaxies (Wise et al. 2008; Greif 
et al. 2008), the initial evolution of H 11 regions therein may 
be similarly confined. This possibility notwithstanding, we 
expect that the spatial resolution that we do achieve suffices 
to track changes in the luminosity emitted in recombination 
lines and in the escape fraction of ionizing radiation, which 
we discuss in the remainder of this Section. 

3.2 Star formation rate indicators 

The luminosity emitted in recombination lines, such as Ha, 
has been found to scale remarkably well with the SFR 
of galaxies at low redshift (e.g. Kennicutt 1983; but see 
Pflamm-Ahenburg et al. 2007). The SFR obtained using 
such relations relies on some knowledge of the IMF of the 
stars which are forming, as well as on the escape fraction of 
ionizing radiation. Fig. 4 shows our calculations of the es- 
cape fraction /cscHii of H l-ionizing photons for each of our 
choices of IMF and total mass in stars, and Fig. 5 shows 
the escape fraction /cscHcIii of He ll-ionizing photons. The 
corresponding luminosities emitted in Hq, Lya, and He 11 
A1640 are presented in Fig. 6. 

As shown in Fig. 4, there is a clear trend toward higher 
H l-ionizing photon escape fractions for more massive stellar 
clusters, with the majority of the ionizing photons escaping 
from clusters with the larger total stellar mass. The escape 
fraction is not, however, independent of IMF; for a given 
total mass in stars, the escape fraction can differ by a wide 
margin. Both the variability in the escape fraction with time 
and the range of values that we find are in rough agreement 
with other recent calculations of the escape fraction of ion- 
izing photons from dwarf galaxies at z 10 (Wise & Gen 
2009; Razoumov & Sommer-Larsen 2009). For a recent cal- 
culation of the escape fraction of ionizing photons from more 
massive galaxies, see Gnedin et al. (2008). 

The breakout of the H 11 region generated by the less 
massive 100 M© star cluster occurs after ~ 1 Myr, lead- 
ing to an escape fraction 0.5 after 2 Myr. In contrast, the 
H II region of the equally massive 25 M© star cluster remains 
confined to the host halo for > 3 Myr, contributing no ioniz- 
ing photons to the IGM. The progress of the initial stages of 
hydrogen reionization, likely driven by star formation in the 
first galaxies (e.g. Loeb 2008) , may thus depend on whether 
these galaxies hosted massive (<; 10 M©) or very massive 
stars (^ 100 Mo)(see also Choudhury & Ferrara 2007). 

The evolution of the H l-ionizing photon escape fraction 
is reflected in the evolution of the luminosity of hydrogen re- 
combination lines, as shown in Fig. 6. Gomparing the panels 
on the left to those on the right, the luminosity in the Lya 
and Hof lines, while generally higher for larger total mass 
in stars, does not scale with the total mass in stars. In- 
deed, owing to the increase in the escape fraction of ionizing 
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Figure 4. The escape fraction of hydrogen-ionizing photons, 
/eacHlli from the host galaxy, each line corresponding to a dif- 
ferent choice of IMF and total mass in stars, as labeled. Note the 
tight anticorrelation between the escape fraction plotted here and 
the luminosity in the hydrogen recombination lines shown in Fig. 
6, demonstrating that the vast majority of the energy emitted in 
hydrogen recombination lines emanates from the dense ionized 
gas within the host halo, as is shown in detail in Fig. 10. 



photons, after ~ 1 Myr the luminosity in hydrogen recombi- 
nation lines from the clusters with greater total stellar mass 
drops below that of the clusters with lower total stellar mass, 
for a given IMF. Overall, because of the temporal evolution 
of the luminosity in a given line, there is no one-to-one rela- 
tionship between the total mass in stars and the luminosity 
in a given recombination line. There is thus likely to be a 
relatively weak correlation between the SFR and the lumi- 
nosity in the hydrogen recombination lines emitted from the 
first dwarf galaxies, owing to the dynamical evolution of the 
photoionized gas and the escape of ionizing radiation into 
the IGM. 

Similar to the case of the hydrogen recombination lines, 
the luminosity in the He ll A1640 line is anticorrelated with 
the escape fraction of He ll-ionizing photons, shown in Fig. 
5. However, different from the case of the hydrogen lines, 
the luminosity emitted in He ll A1640 line is generally much 
more strongly correlated with the total mass in stars, for a 
given IMF. This is due to the low escape fraction of He ll- 
ionizing photons, which is essentially zero for every case 
studied here, except for the case of the more massive cluster 
of 100 Mq stars. With such a high fraction of He ll-ionizing 
photons being balanced by recombinations within the host 
halo, there is a near linear relation between the total mass in 
stars and the luminosity emitted in He ll A1640, making this 
line a potentially much more reliable SFR indicator than hy- 
drogen lines such as Ha. There are slight departures from 
linearity due to the temperature dependence of the emission 
coefficient jAi640, which varies by a factor of ~ 2 over the 
temperature range of the ionized gas in our simulations and 
is generally lower for the hotter H ll regions generated by the 
more massive stellar clusters (Osterbrock & Ferland 2006). 
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Figure 5. The escape fraction of He ll-ionizing photons, /cscHollli 
from the host galaxy, each line labeled as in Fig. 4. For most 
cases, the negligible escape fraction leads to a tight correlation 
between the luminosity emitted in the He II A1640 line and the 
total mass contained in stars, in contrast to the weaker correlation 
for hydrogen recombination lines, as discussed in Section 3.2. 



The luminosity emitted from a galaxy in recombination lines 
depends not only on the stellar IMF, but also on the density 
field of the galaxy and the escape fraction of ionizing pho- 
tons. Therefore, the utility of recombination line strengths 
as IMF indicators hinges on an understanding of the dy- 
namical evolution of the photoionized gas, especially for the 
case of starbursts in the first dwarf galaxies, in which such 
dynamical effects can be most pronounced. 

For the starbursts that we simulate here, the luminos- 
ity of the He ll A1640 emission line relative to the hydrogen 
recombination lines can be read from Fig. 6, while the equiv- 
alent widths of these lines are presented in Fig. 7. Comparing 
the top panels of Fig. 6 to the bottom panels, it is evident 
that the ratio of the luminosity emitted in He ll A1640 to 
that in Ha (or Lya) can be very different depending on the 
IMF. Fig. 7 shows that there is a similar distinction in the 
ratios of the EWs. For the 100 Mq star clusters the luminos- 
ity in He ll A1640 is comparable to that in Ha, while for the 
25 Mq star clusters the luminosity in He ll A1640 is up to an 
order of magnitude lower than that in Hq. However, as the 
escape fraction of H l-ionizing photons increases with time 
for the more massive 25 M© stellar cluster, the luminosities 
in these two lines become comparable, revealing that there 
is some ambiguity in the use of this ratio of line luminosities 
as an indicator of the IMF of young 3 Myr) stellar clus- 
ters. Thus, in some cases dynamical effects may compromise 
the use of this line ratio in distinguishing between clusters 
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Figure 6. The luminosity of the galaxy, as a function of the 
time from the prompt formation of the cluster, in three recom- 
bination lines: Lya {dot-dashed blue), Ha (solid red), and He II 
A1640 (dashed black). The four panels correspond to our four dif- 
ferent choices of IMF and total mass in stars; these are, clockwise 
from top-left: twenty-five 100 stars, two hundred fifty 100 
M(7) stars, one thousand 25 stars, and one hundred 25M 



stars. The luminosities generally decrease with time, as the pho- 
toheating acts to decrease the density of the ionized gas, lowering 
the recombination rate. Note the different evolution of the He ll 
A1640 luminosity as compared to that of the hydrogen recombi- 
nation lines, owing to the lower escape fraction of He ll-ionizing 
photons (see Fig. 5). 



of Pop III.l and Pop III. 2 stars, with typical masses of order 
100 Mq and 10 Mq, respectively. 

The ratio of the observed fluxes in He ll A1640 and 
Ha, as calculated using equation (10), is displayed in Fig. 
8. In this Figure, it is clear that this line ratio is sensitive 
to the IMF, although it is not a constant for each cluster. 
Instead, for clusters in which the escape fraction of H I- 
ionizing photons increases with time dramatically, while the 
escape fraction of He Il-ionizing photons remains roughly 
constant, this line ratio varies with the flux observed in Ha. 
While the ratio of the fluxes is a somewhat ambiguous IMF 
indicator, the clusters with the more top-heavy IMF do con- 
sistently exhibit larger ratios of He ll A1640 to Ha. Nagao 
et al. (2005) present a search for He ll A1640 emission from 
a strong Lya emitter at z = 6.33, flnding an upper limit 
for the ratio of He ll A1640 to Lya. Assuming a standard 
value of 0.07 for the ratio of luminosity in Ha to that in Lya 
(Osterbrock & Ferland 2006), we show in Fig. 8 the upper 
limit that these authors report (see also Dawson et al. 2004) . 
Although a weak upper limit, it is clear that observations 
with only slightly greater sensitivity will allow to differenti- 
ate between the flux ratios predicted here for massive and 
very massive Pop HI IMFs. 

Comparing the EW of Ha in the four panels of Fig. 7, 
it is clear that it is not strongly dependent on the IMF or 
on the total mass in stars, varying by at most a factor of 
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Figure 7. The rest frame equivalent widths, = W / (I-I-2;), 
where W is the observed EW, as a function of time, of the same 
three recombination lines shown in Fig. 6. For comparison, in each 
panel we plot the observed EWs of galaxies from two different 
surveys: the dotted line at 400 A denotes the median EW of Lya 
emitters detected at 2; = 4.5 in the LALA survey (Malhotra & 
Rhoads 2002), while the dotted line at 145 A denotes the average 
EW of the six Lya emitters at z > 6 detected in the Subaru 
deep field (Nagao et al. 2007). Note that the Lya EWs that we 
compute are upper limits, as scattering in a neutral IGM has not 
been accounted for. 



three between each of the cases. While showing more vari- 
ation between the four cases, the EW of Lya also shows 
considerable ambiguity as an IMF indicator, its maximum 
value varying by about a factor of three between each of 
the four cases. This insensitivity of the Lya EW to the IMF 
arises from two effects. Firstly, the stellar continuum lumi- 
nosity La,* increases in a similar manner as the number of 
ionizing photons from the massive (25 Mq) IMF to the very 
massive IMF (100 M©). This acts to keep the EW, roughly 
the ratio of the two, relatively constant. Secondly, while the 
luminosity in Lya decreases with the increasing escape frac- 
tion of ionizing photons for the more massive clusters, the 
continuum luminosity remains largely unchanged, leading to 
a decrease in the EW with time for these clusters. We note 
that the Lya EWs presented here are only upper limits, as 
we have not accounted for scattering of Lya photons in the 
IGM (see e.g. Dijkstra et al. 2007). 

The EW of He 11 A1640 is a more definitive indicator of 
IMF, being higher for the clusters of 100 Mq stars than for 
the clusters of 25 Mq stars, regardless of the total mass in 
stars or of the age of the cluster (up to J; 3 Myr). As with 
the utility of He 11 A1640 as a SFR indicator, this largely fol- 
lows from the generally low escape fraction of He ll-ionizing 
photons. 

For comparison with observed galaxies, we plot in Fig. 7 
the two observational results: the median Lya EW of galax- 
ies detected in the Large Area Lyman Alpha (LALA) survey, 
W° ~ 400A, and the average EW of six galaxies observed 
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Figure 8. The ratio of the integrated fluxes in He II A1640 and 
Hcf, F;^ig4o/FHc«i as a function of time, for each of the four clus- 
ters simulated here, as labeled. The dashed horizontal line de- 
notes the upper limit of this ratio for the strong Lyo emitter 
SDF J132440.6-I-273607 at z = 6.33, as reported by Nagao et al. 
(2005). Similar upper limits for Lya emitters at z = 4.5 have been 
reported by Dawson et al. 2004. 



Figure 9. The observed fluxes, /, as a function of time, of the 
same three recombination lines shown in Fig. 6, for the galaxy we 
simulate at z ~ 12.5. The fluxes are normalized to what would 
be observed with a spectroscopic resolution of .R = 1000, and are 
computed using equation (10) assuming that the galaxy appears 
as an unresolved point source. Note that the flux in Lya is an 
upper limit, as the present calculation does not take into account 
scattering in a neutral IGM. 



at > 6 in the Subaru deep field, ~ 145A (Nagao et al. 
2007). The large LALA EWs are comparable to what we find 
for primordial dwarf galaxies, although the LALA galaxies 
likely do not host Pop III star formation (but see Jimenez & 
Haiman 2006). The detection of an EW of the He ll A1640 
line lOA would be a stronger indication of a galaxy host- 
ing primordial star formation, as shown in Fig. 7, although 
none has been found as of yet. We note that observed Ly- 
man break galaxies at z ~ 3 have been found to have He ll 
A1640 EWs of ~ 2A (Shapley et al. 2003), consistent with 
what is expected for Wolf-Rayet stars formed in starbursts 
(see Schaerer & Vacca 1998; Brinchmann et al. 2008). 

3.4 Detectability of Recombination Radiation 

In Fig. 9, we present our predictions for the observable re- 
combination line fluxes, for each of the stellar clusters that 
we simulate. Fig. 10 shows the surface brightness in Ha as 
observed on the sky, for the two more massive stellar clus- 
ters, which each have a total mass in stars of 2.5 x 10* M©. 
The fluxes in each plot, largely determined by our choices 
for the total stellar mass, are calculated using equation (10). 
While the larger H ll region generated by the more massive 
stars encompasses more dense gas, creating more widely dis- 
tributed emission in Ha, as shown in Fig. 10, the highest 
flux per square arcsecond is in the central region of the halo 
hosting the less massive stars. This is due again to the less 
dramatic dynamical response of the gas to photoheating, 
leading to higher densities, and thus higher recombination 
rates. Due to this effect, the highest fluxes are generated 
just after the birth of a stellar cluster, as shown in Figure 



9, when the density of the photoionized gas is still high, not 
having had time to expand in response to the concomitant 
heating. Indeed, Fig. 9 shows that the flux in Ha from the 
more massive 25 Mq star cluster may reach ^ 10 x (7?/1000) 
njy before the breakout of the H ll region. Catching the flrst 
galaxies when still in the earliest stages of their initial star- 
bursts, within the first few 10^ yr, is thus likely to provide 
one of the best chances for observations of purely primordial 
stellar populations in the early Universe. 

Aboard the JWST, the Ha line would be observed with 
the Mid-Infrared Instrument (MIRI). Its pixel size of 0.1 
arcsec would not resolve the brightest portions of the galax- 
ies that we simulate, which, as shown in Fig. 10, are of order 
0.01 arcsec. With a resolution capability of .R = 3000 the 
MIRI has a sensitivity of <; 200 nJy for a signal-to-noise of 

10, in exposures of ~ lO'' s (see Panagia 2004), making it un- 
able to detect even the brightest galaxies that we simulate, 
the flux in Ha of these being 20 nJy for R = 3000. 

With a greater sensitivity of ^ 100 nJy (Panagia 2004), 
the Near Infrared Spectrograph (NIRSpec) operates in the 
wavelength range 0.7 to 5 /im, allowing it to possibly detect 
Lya out to 2: ~ 40 and He 11 A1640 out to 2 ~ 30. However, 
the flux in He 11 A1640 is always below that in Ha and, hence, 
is too low to be detected. The Lya line, with the highest flux 
of the three recombination lines shown in Fig. 9, is also not 
directly detectable, with a flux falling well below the ~ 100 
nJy sensitivity limit of NIRSpec. Furthermore, although the 
luminosity in Lya is always intrinsically higher than that in 
Ha, before reionization the observable flux in Lya may be 
dramatically decreased due to scattering in the neutral IGM 
(e.g. Dijkstra et al. 2007) . Because we do not account for this 
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Figure 10. The flux in Ho, /hqi per square arcsecond, emitted from a primordial dwarf galaxy, as observed on the sky at 2: = 0, 
assuming a spectroscopic resolution of R = 1000. Shown here are the two most massive of the four stellar clusters that we simulate, 
one containing 25 Mq stars (bottom panels), the other containing 100 Mq stars {top panels). From left to right, the clusters are shown 
at 10^ yr, 1 Myr, and 3 Myr after formation. Note that the emission is concentrated in the densest photoionized regions, the filaments 
around the galaxy (compare the top panels to Fig. 2) and especially the dense gas within the inner ~ 100 pc of the galaxy. The highest 
total fluxes occur at the earliest times, before the H n region has broken out; hence, the youngest stellar clusters are the most readily 
observed. 



effect in the present calculations, the Lya fluxes presented 
here are only upper limits. 

The Near Infrared Camera (NIRCam) aboard the 
JWST, which will be used to conduct deep surveys designed 
to detect the first galaxies, will be capable of detecting point 
source fluxes as low as ~ 3.5 nJy at a signal-to-noise of 10 
for a 10® second exposure (e.g. Gardner et al. 2006). With 
a resolution of 0.03 arcsec per pixel, the NIRCam would 
also not quite resolve the galaxies that we simulate. We can 
evaluate the possibility that NIRCam may detect them as 
point sources, however, by estimating the continuum flux of 
the galaxies as observed at ~ 2 /xm. As can be read from 
Figs. 6 and 7, the continuum flux, oc Lai640 / W^ai640, varies 
only by a factor of a few between Lya and He 11 A1640. Thus, 
for simplicity we assume that the continuum is roughly flat 
and calculate the specific continuum flux, as observed at z 
= 0, as 

r -^A1640 ^A164(l(l + ^) 

~ 0.03nJy , (16) 

where VK^m4Q is the equivalent width in the rest frame of the 
galaxy, as defined in equation (9) . This flux is well below the 
sensitivity limit of the NIRCam, and so we conclude that 
detection of the continuum radiation from the galaxies we 
simulate would also be undetectable. We note, however, that 



under favorable circumstances, gravitationally lensed emis- 
sion from a primordial galaxy undergoing an only slightly 
more luminous starburst may be detectable with the JWST, 
given that lensing can boost the flux by a factor of order 10 
(e.g. Refsdal 1964; Stark et al. 2007). 

The first dwarf galaxies could be more luminous than 
we find here if the efficiency of star formation esF, defined 
as the fraction of the total baryonic mass in the galaxy con- 
tained in stars, is larger than what we have assumed in our 
simulations, where our choices for the total mass in stars cor- 
respond to modest values of esF ~ 10"^ -10"". A larger efii- 
ciency (esF ~ 10"^), for the top-heavy IMFs considered here, 
would yield a cluster observable by NIRCam (e.g. Gardner 
et al. 2006). However, as we have demonstrated, the much 
higher ionizing flux from a IQp M© cluster of massive 
primordial stars would induce a strong hydrodynamic re- 
sponse which would lead to a rapid decline in the luminos- 
ity emitted in recombination radiation. Thus, even if such 
clusters can be identified by their continuum emission, the 
detection of recombination radiation, and with it informa- 
tion about the stellar IMF, may be beyond the capabilities 
of the JWST. Furthermore, the formation of such a massive 
cluster of primordial stars may face an impediment due to 
the strong radiative feedback within the cluster itself. Re- 
cent simulations of clustered star formation in the present- 
day Universe suggest that radiative feedback infiuences the 
fragmentation behavior of the gas and possibly lowers the 
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overall star formation efficiency (see Krumliolz et al. 2007; 
Bate 2009). However, the situation is by no means clear (e.g. 
Dale et al. 2005, 2007). 

It is possible that more massive (10^ - 10^° M©) pri- 
mordial galaxies form at 2 ^ 12, or form at lower redshift, 
making their detection feasible. However, such more mas- 
sive, and therefore more luminous, galax;ies are likely to also 
be more chemically evolved, and so may already be domi- 
nated by Pop H star formation. Thus, it may be that the 
galaxies which host pure Pop HI starbursts, such as those 
we study here, will remain out of reach of even the JWST, 
although this critically depends on the poorly constrained 
process of metal enrichment in the early Universe (see e.g. 
Pan & Scalo 2007; Tornatore et al. 2007; Cen & Riquelme 
2008; Johnson et al. 2008). We emphasize, however, that the 
dynamical effects studied here are likely to play a role even 
in the more luminous galaxies that will be detected, and are 
important to take into account in evaluating observations 
meant to constrain the SFR or the IMF. 



4 SUMMARY AND CONCLUSIONS 

We have presented calculations of the properties of the 
recombination radiation emitted from a primordial dwaxf 
galaxy at z ~ 12.5, during the initial stages of a starburst. 
Our cosmological radiation-hydrodynamical simulations al- 
low us to track the detailed dynamical evolution of the emit- 
ting gas in the central regions of the galaxy, and thus to 
study its effect on the emerging radiation. The goal of this 
study has been to determine the observable signatures of the 
initial starbursts in the first galaxies. In particular, we have 
aimed to find reliable indicators of the star formation rate 
and of the stellar IMF. 

Owing to the escape of H i-ionizing photons into the 
IGM, we find only a weak correlation between the total mass 
in stars and the luminosity in hydrogen recombination lines. 
This suggests that Lya and Ha, despite the high luminosity 
in these lines, may not serve as strong indicators of the SFR, 
unlike in the low-redshift Universe (e.g. Kennicutt 1983). 
The He ii A1640 line may be a more effective SFR indicator, 
as the luminosity in this line scales more closely with the 
total mass in stars, due to the lower escape fraction of He ii- 
ionizing photons. 

We confirm that the ratio of He ll A1640 to either Lya 
or Ha can be used as an IMF indicator, although its utility is 
compromised in some cases by the unequal escape fractions 
of H I- and He ii-ionizing photons. The most robust IMF 
indicator, in terms of distinguishing between populations of 
massive (^ 10 M©) and very massive (S; 100 M©) Pop HI 
stars, is the EW of He ll A1640, as it is consistently higher 
for the more massive stars regardless of the total mass in 
stars. We note that while in principle the radiation emitted 
by a central accreting black hole (BH) in a primordial dwarf 
galaxy could introduce complications for using He ll A1640 
as IMF indicator (e.g. Tumlinson et al. 2001), recent work 
suggests that BH accretion is inefficient in the early Universe 
(e.g. Johnson et al. 2007; Pelupessy et al. 2007; Alvarez et 
al. 2008; Milosavljevic et al. 2008, 2009). Such miniquasar 
activity may thus not result in appreciable observable radi- 
ation for the first ^ 10* yr. 

In cases where the gas in primordial dwarf galaxies is 



unable to cool efficiently due to the photodissociation of 
H2, it has been argued that, instead of a stellar cluster, a 
BH with mass ^ 10* M© may form by direct collapse (see 
Bromm & Loeb 2003; Lodato & Natarajan 2006; Begelman 
et al. 2006; Spaans & Silk 2006; Regan & Haehnelt 2008). 
The fraction of dwarf galaxies that host the formation of 
such BHs is likely quite small, since the photodissociating 
flux required to suppress fragmentation is very high, well 
above the cosmological average at 2; ^ 10 (Dijkstra et al. 
2008). In the present work, we have chosen to focus on the 
presumably more common case that a stellar cluster forms 
(but see Begelman & Shlosman 2009). We do note, however, 
that stellar mergers within dense clusters, such as those we 
consider here, may lead to the formation of ~ 10^ M© BHs 
(Omukai et al. 2008; Devecchi & Volonteri 2009). 

In terms of the detectability of the recombination ra- 
diation from the first galaxies, we have shown that due to 
the dynamical response of the gas to photoheating, a top- 
heavy IMF or a high star formation efficiency can be self- 
defeating, leading to a decrease in the line luminosity of the 
galaxy. We conclude that the detection of purely primordial 
dwarf galaxies at « ^ 10 is likely to be beyond the capa- 
bilities of the JWST, although their detection may be just 
possible if the galaxies are strongly lensed. More luminous, 
1O^-1O'"M0 (total mass), galaxies may thus be detected by 
the JWST; however, such more massive galaxies are likely to 
be already chemically enriched due to previous episodes of 
star formation in their progenitors. These systems are then 
expected to host Pop II star formation, or a composite of 
Pop II and Pop HI, depending on the still poorly understood 
mixing of heavy elements in the first galaxies (e.g. Karlsson 
et al. 2008). 

Our results demonstrate how the radiation emitted from 
the first galaxies depends on the hydrodynamic effects of the 
photoionization from clusters of massive stars. This is com- 
plementary to the results of previous studies (e.g. Schaerer 
2002, 2003), which highlight the evolution of the emitted ra- 
diation owing to the aging of a stellar population. Clearly, 
both effects must be considered in future work. 

The initial staxbursts of the first galaxies may constitute 
the formation sites of the only metal-free stellar clusters in 
the Universe, since after the first several Myr supernova feed- 
back can quickly enrich the galaxy with metals (e.g. Mori et 
al. 2002; Kitayama & Yoshida 2005; but see Tornatore et al. 
2007, Cen & Riquelme 2008). Also, a large fraction of the 
first dwarf galaxies, with masses of order 10* Mq, may al- 
ready form from metal-enriched gas (Johnson et al. 2008; see 
also Omukai et al. 2008); it is an important open question 
what fraction of dwarf galaxies forming at « ^ 10 are pri- 
mordial when they begin forming stars. Future observations 
of those first dwarf galaxies that do host primordial star 
formation offer one of the few opportunities for constraining 
the primordial IMF. 
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